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Abstract 


Initial boundary value problems describing in-situ leaching of rare metals (ura- 
nium, nickel, etc.) with an acid solution are considered. We assume that the 
solid skeleton is an absolutely rigid body. Next we describe the physical process 
in the pore space (that is at the microscopic level, where a characteristic size is 
about 5-20 microns) by the Stokes equations for an incompressible fluid coupled 
with diffusion-convection equations for acid concentrations and chemical reaction 
products. Since the solid skeleton changes its geometry during dissolution, the 
boundary ”pore space - solid skeleton” is unknown (free) and the pore space has 
a special periodic structure. To derive a macroscopic mathematical model (char- 
acteristic size meters or tens of meters) for this special periodic structure we use 
Nguenseng’s two-scale convergence method.As a result we prove the existence and 
uniqueness theorem for the obtained macroscopic mathematical model. 
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1 Introduction. 


In the proposed manuscript we rigorously derive a macroscopic mathematical 
model of in — situ leaching of rare metals on the base of homogenization in 
structures with a special periodicity and prove theorems of the existence and 
uniqueness of the classical solution to the initial - boundary value problem 
for the corresponding system of differential equations. 

The extraction of rare metals by leaching is very important task of the 
national economy. Natural deposits of uranium, nickel and other rare metals 


*Moscow State University of Civil Engineering, Moscow, Russia 
e-mail: anvarbek@list.ru. 


Figure 1: General scheme of the in-situ leaching process 


are complex geologically heterogeneous objects. Inhomogeneity means that 
the properties of an object of interest change in space. Analyses of wells and 
cores show that the geological properties (porosity, permeability, etc.) of ore 
bodies are heterogeneous even within a single deposit. Very often insufficient 
taking into account the consequences of inhomogeneities at the stage of op- 
eration planning becomes obvious too late, when the acid solution uploaded 
to soil through injection wells appeared far from the intended location. In 
addition, an important role is played by the concentration of the injected 
acid, the injection modes of acid solutions and other factors. 

Hence, understanding the dynamics of fluids in heterogeneous porous me- 
dia and the mechanism of dissolution of rocks by acids is of fundamental im- 
portance for the effective management of rare metal mining. This is achieved 
by creating a prototype of a hydrodynamic simulator for the ore body based 
on the corresponding mathematical models, which allows to optimize the en- 
tire technological process. A hydrodynamic simulator of an ore body is a 
complex, consisting of a scale of mathematical models of a physical process 
(prototype of the hydrodynamic simulator), digital characteristics of the ge- 
ometric and physical properties of the solid skeleton, and a set of computer 
programs, that allow one to visualize the physical process and determine the 
dynamics changes in the main characteristics of the mathematical model. 

Currently, a large range of mathematical models exist to describe the 
dynamics of rock leaching, which may describe the physical processes un- 
der consideration, but only at the Macroscopic Level (sce [1] - [3] and 
references there). Unlike Microscopic Models (the characteristic size is 
approximately tens of microns) in Macroscopic Models, the characteristic 
size is meters or tens of meters. Because of this, these models do not dis- 
tinguish between the microstructure of a continuous medium, since in such 
a model at each point the medium contains both solid skeleton and liquid in 
the pores or cracks of this skeleton. 


All such models are built on the same principle. Fluid dynamics is usually 
controlled by the Darcy’s system of filtration, or some modification of it. The 
equations describing the migration of acid and chemical reaction products 
simply postulated and, roughly speaking, are some modifications of diffusion- 
convection equations for the corresponding concentrations. The main thing 
in these postulates is the type of coefficients of the equations. 

Exactly here there is a great variety of models, depending on the tastes 
and preferences of its authors. It is quite understandable, since the main 
mechanism of the physical process is focused on the unknown (free) boundary 
between the pore space and the solid skeleton and not spelled out in any 
way in the proposed macroscopic models. This is where the rocks dissolve, 
changing the concentration of the injected acid, and this is where products 
appears inside the carrier liquid. Moreover, during the process, the geometry 
of the pore space (geometry of the boundary separating solid skeleton and 
pore space) changes in time and space. 

These fundamentally important changes occur at the microscopic level, 
corresponding to the average size of pores or cracks in rocks, while all of 
proposed macroscopic models operate with completely different orders of 
scales (tens of centimeters or meters) and, therefore, do not distinguish free 
boundaries, nor features of the interaction of acid with rocks, which explains 
the wide variety of macroscopic mathematical models. 

The authors of such models simply do not have an exact method for 
describing physical processes at the microscopic level based on the funda- 
mental laws of the classical continuum mechanics and chemistry, nor the 
possibility to take into account the microstructure in their macroscopic mod- 
els. Therefore, they have to restrict themselves with some certain speculative 
considerations. 

Due to this, a natural question arises. If there are several macroscopic 
models describing the same physical process under the same conditions, 
which of them most adequately reflects this process? Where is the criterion 
of the adequacy here? It doesn’t make sense to talk about an experiment, 
since each such model has enough free parameters that are not related to 
reservoir geometry (e.g. porosity), or to the physical characteristics of the 
process (such as the viscosity of the filtered liquids, or the physical proper- 
ties of the solid skeleton and the like). So, with the variation of these free 
parameters one may get a match with any experiment. 

R. Burridge, J. B. Keller [4] and E. Sanchez-Palencia [5] were the first 
who explained that the exact description a filtration of liquids and seismic 
in rocks at the macroscopic level is possible if and only if: 

(a) the physical process under consideration is described at the micro- 
scopic level by equations of Newtonian classical continuum mechanics (exact 
model); 

(b) a set of small dimensionless parameters is selected; 


(c) the macroscopic mathematical models is an exact asymptotic limits 
(homogenization) of exact mathematical models at the microscopic level, 
when the selected small parameters tend to zero. 

Various special cases of exact macroscopic models of acoustics and fluid 
filtration in rocks have been investigated by many authors: Gilbert-Lin [6], 
Ferrin, Mikelic [7], Levi [8], Sanchez-Hubert [9], Zhikov, Kozlov, Oleinik [10], 
Zhikov [11], Pastukhova [12], Bachvalov, Panasenko [13] and others. All 
these authors used various methods of homogenization and, as a rule, the 
resolution to each of them was a difficult task required considerable effort 
and ingenuity. 

Everything changed after appearance of G. Nguetseng’s paper[14], where 
the author proposed the Method of Two-Scale Convergence in peri- 
odic structures. What used to be an art has become an ordinary routine, a 
reference to the method. So the homogenization theory has ceased to be an 
independent part of the mathematical analysis (or the theory of differential 
equations) and main efforts in homogenization have moved from theory to 
applications in mechanics, physics, biology, etc. 

Since the problems in periodic structures have been well studied, the 
attention of researchers began to attract more complex problems in Me- 
dia with special periodicity, set by characteristic functions of the form 
x° (x,t) = x(a,t, —), where x(a, t, y) — 1 — periodic by variable y function. 

The most fully exact models of physical processes in poro-elastic me- 
dia were investigated by A. Meirmanov [15]. He rewrote these models in a 
special dimensionless form and fixed a certain Set of criteria (dimension- 
less coefficients of differential equations) that are responsible for the type 
of physical process (filtration, acoustics, hydraulic shock, etc.). Meirmanov 
classified accurate mathematical models at the microscopic level and using 
the Nguetseng multi-scale convergence method obtained exact asymptotic 
limits of these models that adequately describe the physical processes under 
consideration at the macroscopic level. 


In our manuscript, the value « = —, is taken as a small dimensionless 


parameter. Here | is characteristic pore size and L is the characteristic size 
of the physical domain under consideration. The dimensionless parameter 
a’, characterizes the viscosity of the filtered liquid: 


a= a 
ML gt po’ 


where 7 is the characteristic duration time of the physical process, po is 
the density of water, g is the acceleration of gravity and yz is the dynamic 
viscosity of the liquid. 

As we have already noted, the derivation of macroscopic mathematical 
models should be based on the most accurate mathematical model of the 


physical process at the microscopic level, described by the laws of classi- 
cal continuum mechanics. In our description, these are the Stokes equa- 
tions in the pore space (2; for an incompressible viscous fluid with dynamic 
characteristics v*(a,t)( Velocity vector (in the liquid component), p* (a, t) 
(Pressure (in the liquid component) and the diffusion-convection equations 
for Concentration acid c*(x,t) and Concentrations of chemical re- 
action products c;(x,t), 7 =1,...,n. 

The solid skeleton Q£ is assumed to be Absolutely rigid in which the 
velocities of the medium are zero. 

Differential equations are supplemented by boundary conditions at the 
given boundaries, initial conditions and strong discontinuity conditions on 
the unknown (Free) boundary I* ”solid skeleton-pore space” in the domain 
Q = 05 UT* US [16],arising from the conservation laws of classical mechanics 
and an additional condition on the free boundary 


DSO C50 = const. > 0, (1.1) 


postulated in theoretical chemistry and allows us to find a free boundary [17]. 

Here D%, is the velocity of the boundary I in the direction of the unit 
normal N* to the boundary I*, outward to the domain 0%, a is a given 
constant. 

The formulation of the problem was previously given in the monograph 
[18], but there were no exact results about the existence of any solution in 
this monograph. 

Note that the physical process we are considering has a rather long dura- 
tion (the filtration rate of the liquid is several meters per year). Therefore, 
the most interesting are the theorems of the existence of solutions of the cor- 
responding initial - boundary value problems globally in time. On the other 
hand, due to the strong nonlinearity of problems with free boundaries [19], 
it is usually not possible to prove any result globally in time for mathemat- 
ical models at the microscopic level. That is, possible results can only be 
theorems of the existence of a generalized or classical solution of the initial - 
boundary value problem for a system of differential equations describing the 
leaching process at the macroscopic level. 

But how can we obtain macroscopic mathematical models if we know 
nothing about the existence of solutions to microscopic mathematical models, 
the limit of which they should be? 

To get around these difficulties, we will use the fixed point theorems [20]. 
To do this, the structure of the pore space is fixed, given by the character- 
istic function x(a,t, y). As we have already noted, in the case of a general 
provision,to solve the emerging problem is almost impossible. Therefore, it 
is reasonable to limit ourselves to the simplest cases. For example, when a 
non-negative function r(x, t) from some set Mtr defines a characteristic func- 
tion of the pore space y(r, y). Next, for a fixed r € Mr, the initial boundary 
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value problem B¥*(r) is considered in a given domain 4(r) for determining 
the main characteristics of the medium (velocity vector, pressure and acid 
concentration), which is the initial boundary value problem A* without an 
additional boundary condition (1.1). In this auxiliary problem, for a fixed 
€ > 0, the solid skeleton is a union of disjoint sets sufficiently close to balls 
of radius er, slowly decreasing in volume, which simplifies the geometry of 
the original pore space and allows us to prove the existence of approximate 
solutions. 

To understand what should be theHomogenized Problem H(r) of the 
problems B*(r), a formal homogenization of the problem A* is performed 
beforehand. The conditions for the existence of such homogenization are 
formulated in the lemma 4.1. 

If r°(a,t) defines the structure of the solid skeleton (pore space) in the 
problem A® and rz > r* at ¢ + 0, then homogenize problem H(r*) of the 
problem B*(r*) = A® should coincide with the homogenization H of the 
problems A* without homogenization of the boundary condition (1.1). It is 
clear that the homogenization of the boundary condition (1.1) with a given 
pore space structure r(a,t) forms a problem operator whose only fixed point 
r*(a,t) determines the required unique homogenization H of the problem A*. 

To solve the problem H(r) we, first of all, have to solve the problem Bé(r) 
and then find its homogenization H(r) for e« — 0. 

The linear problem B*(r)decomposes into a sequential solution Dynamic 
problem B*(r) on the definition of dynamic characteristics v°(a,t) and 
p*(a,t) and Diffusion-convection problems B*(r) for finding the acid 
concentration c*(a,t) and concentrations of products of chemical reactions 
cf (a, t),j=1,..,n. 

Due to the linearity of these problems, the existence and uniqueness of 
a weak solution to each of them follows from the corresponding a priori 
estimates and known methods for solving linear differential equations. For 
example, the Galerkin’s method [21]. 

The next step is to homogenization the resulting mathematical model 
B*(r). To homogenize this problem, we use the Nguetseng’s two-scale con- 
vergence method [14], which makes it quite easy to get a mathematical model 
H(r). But since this method was developed only for averaging functionals, 
we needed to write down the original mathematical model in the form of a 
system of integral identities equivalent to the original system of differential 
equations. 

If the integral identities equivalent to the dynamic Stokes equations, as 
well as the integral identities for the diffusion-convection equation with clas- 
sical boundary conditions, are well known, then a similar identity for acid 
concentration, including the diffusion-convection equation and boundary con- 
ditions on the free boundary, has been unknown until now. 

Equivalent expression of differential equations in the form of integral iden- 


tities is a common and rather difficult task for free boundary problems. The 
happy exception was the Stefan problem [22], [23], describing phase transi- 
tions in pure (without impurities) media. Such, for example, as ” water — ice” 
or chemically pure metals. The authors of the cited works O. A. Oleinik and 
S. L. Kamenomostskaya (S. Kamin) managed to reformulate the problem so 
that the equivalent formulation in the form of an integral identity, in the case 
of the existence of a classical solution, contained both the equation of thermal 
conductivity outside the free boundary and the condition itself on the free 
boundary. This approach allowed them to prove quite simply the existence 
and uniqueness of the weak solution to the integral identity under minimal 
conditions for the smoothness of the solution. At the same time, in the case 
of the existence of a classical solution of the Stefan problem, the latter must 
coincide with the weak solution. The question of the existence of a classical 
solution to the one-phase Stefan problem remained open until 1975 [24]. The 
existence of a classical solution to the two-phase Stefan problem was proved 
in 1979 [19]. 

Considering the above, in the mathematical model of in — situ leaching 
proposed by us, it was first of all very important to find such an equivalent 
formulation of the problem in the form of a system of integral identities that 
would require minimal smoothness of solutions to the problem, which was 
done. 

At the same time, for the fluid moe v* and the acid concentration 


, defined only in pore space 04 7(r -U OF(r), it was necessary to find 


an extension of these solutions from the dorian of definition onto domain 
Qr = Q~x (0,7) while preserving their best differential properties. To do this, 
we used the results on the extension of functions, formulated in Lemmas 2.8 
and 2.9 § 2, so that the extensions v “(a,t), p°(a,t) and ¢°(a, t) of functions 
v (x,t), p°(a,t) and c*(a,t) satisfy the integral identities of conservation 
of momentum, mass and acid concentration, equivalent to the corresponding 
differential equations together with the boundary and initial conditions. 

The moving boundary I*(r) creates the main problem when deriving 
a’priori estimates to solutions to the problem B*(r). Standard methods to 
getting a priori estimates by multiplying a differential equation by solutions 
(or a combination thereof) of these equations in free boundary problems often 
do not lead to the desired result. In our situation, this method helped only 
partially. 

Similarly, the derivation of a priori estimates of solutions to the diffusion 
— convection problem B*(r), as a rule, is based on the principle of the max- 
imum inherent in standard initial - boundary value problems for parabolic 
equations. Unfortunately, there is no maximum principle in the initial for- 
mulation of our problem. Unfortunately, there is no maximum principle in 


the initial formulation of the problem. 

All of the above required the introduction of the Modified diffusion 
— convection problem B¥*(r),which coincides with diffusion — convection 
problem B*(r), if the concentration of acid in modified diffusion — convection 
problem B*(r) satisfies the constraints of the standard maximum principle 
(Osa): 

A’priori estimates of weak solutions (that is, solutions to the correspond- 
ing integral identities), as a rule need a special selection of trial functions 
in the integral identity and integration in parts. For the latter, sufficient 
smoothness of the boundary of the pore space 24(7) (the domain filled with 
liquid) is necessary. The smoothness of the boundary 00%4(r) is defined by 
r € Mp. This simple fact is central one to the derivation of a’priori estimates. 

Let’s briefly focus on the structure of the article. 

In § 2 well-known facts and definitions are given, as well as new results 
on the compactness of sequences in a non-periodic structure (Theorem 2) 
and results on the extension of functions defined in a pore space 24.-(7) with 
a given structure, defined by the function r(a,t) are proved (Lemmas 1.2.8 
and 1.2.9). 

§ 3 is devoted to the formulation of an initial boundary value problem 
describing a physical process at the microscopic level. 

In § 4 the formal homogenization of the problems A® and B*(r) are con- 
sidered. 

To get integral identities equivalent to differential equations supplemented 
with the corresponding boundary and initial conditions we assume the nec- 
essary smoothness of the solutions to the original problem. These integral 
identities allow us to correctly define the weak solutions to problems A* and 
Bé(r). 

In Lemmas 4.1 —4.3 we formulate the necessary conditions for homoge- 
nization the problems A*® and B*(r) and formally derive homogenized models 
H(r) and H or the limiting velocity, pressure and acid concentration. 

In Lemma 4.2 under assumption (4.18) we derive integral identity (4.20), 
equivalent to the boundary condition (1.1), and find the homogenization 
(4.19) of this boundary condition. On the base of equality (4.19) we construct 
the problem operator 


P(r) = 6 7 c(w, T)dr. 


The fixed point of this operator define the desired structure of the pore space 
and auxiliary operators F°(r) = c(a,t), F°(r) = v(a,t) F(r) = p(a,t), 
where {v, p, c} are solutions to the problem H(r). 

The constant @ is defined by equality (4.18). 

In Lemma 4.3 we prove the correct definition of these operators F°(r), 
F’(r) and F(r). 

In § 5 we formulate the main results of the manuscript. 


In § 6 we prove Theorem 3. The crucial moment here is a’priory estimates 
for the problem B*(r). 

In § 7 we prove Theorem 4 on existence of weak and classical solutions 
to the problem H(r). 

§ 8 is devoted to the prove of Theorem 5. In this theorem we analyze 
properties of opderators F°(r), F°(r) and F?(r)and the principal operator 
F(r). 

We show that operator F(r) is a Lipschitz continuous one with the cor- 
responding constant bounded by some linear function of 7’. This fact proves 
the existence of the unique fixed point r*(a,t) locally in time. 

The last statement means the uniqueness of the solutions to the mathe- 
matical model H. 

Finally, using smoothness of the solutions to the problem H(r) we prove 
the correctness of the mathematical model H for any T' > 0. 

In our paper, we used the notation adopted in [25]and in [26]. 


2 Auxiliary statements. 


2.1 Notations. 
2.1.1 Dimensionless parameters. 


We have already defined dimensionless parameters a,, c¢ in the introduction. 
Diffusion of acid and chemical reaction products is characterized by di- 
mensionless diffusion coefficients 


respectively, acids and products of chemical reactions. 
We will assume that a,, depends on the small parameter ¢: a, = aj, and 


there are finite or infinite limits 
é 


lima’ = uo, lim & = py. 
e>0 i Ho, e>0 e2 


The coefficients cr, dp a; are fixed and do not depend on the small parameter 
é. 

As 0, we denote the dimensionless density of the solid skeleton, related 
to the density of water pp and as gf — dimensionless density of the liquid 
component related to the density of water po. 


2.1.2 Domains and boundaries. 


Q Cc R® is a bounded domain with peace piecewise smooth boundary S, 


§ = 00,5 = SUSUS?, Sins® =6,1 =1,2, Sins? =0, Sh = 
Si x (0,T) CR‘, j =0,1,2,3, Ar =Qx (0,T) CR. 
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The boundary 5° C R? is impermeable to liquid in the pore space,, the 
boundary S' C R® simulates injection wells and boundary S? C R? simulates 
production wells. 

For simplicity of presentation, in what follows we assume ( the unit cube, 


1 1 1 
|x2| < gh Ute (22 = £5; |x1| < 9? |x3| < 


. ee “2 Fe “ 7 1 1 il 

af St = 4@ 2% = 2 |x2| < 2 |x3| < 5h 3? = {x 2 fy = ~ 5? |x2| < 
95? |x3| < } 

2 2 


1 
Let r= ret), 0 <r(et) < 5° (a,t) € Qr, be a given function, that 
determines the structure of the pore space Q%(r) C Q and ve oe skeleton 


OS(r) C Qin C and the structure of the pore space 0% 7(r) -U OF (r) C Op 


at 


and the solid skeleton Q5 (r) + OQF(r) C Op in Nr. 
Let also the boundary I*( = 905( nc ke OQ §(r) divides liquid and solid 


components in 2 and ['}.(r is I(r) in Qr. 


2.1.3 Structure of pore space. 


In what follows all functions of the type ®(a,t, y), where (w@,t) EQ y € R® 
are considered 1 - periodic in variable y: 


O(a,t,y) = O(a,t,<(y)), y= [ll] +s(y), llyll = (mld, (yell, ply i 


Here ||a|] is the integer part of the number a. 
For fixed ¢ > 0 the pore space f(r) and the solid skeleton 0{(r) are 
defined by 1 — periodic in variable y characteristic function x(r, y) as 


OF(r) = Int{x €Q: x°(x#,t) = 1}, 
OS(r) = Int{a €Q: x°(a,t) = 0}, OF(0) = 0, 5 (0) 
sen(r(ax,t) — |s(y) 
2 


= i 
— 
~ 
Nw 
nN 


x'(2,t) = x(r(@,0), =), x(ry) = 
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N=|(JO™, OF S{eeO:e=ek fey}, OFF =0** x'(0,7), 
OG) = ae DCPS: ORS: ae OO", 


OF@) =U) (r), OFF) = Yat 
If r(ek,t) =0, ORG (ean Mee ca Chk Oh 


for all k = (ki, ko,k3), ki, k2,k3 € Z (integer numbers) and for all y = 
(yi, Y2, Y3) ea l 
Here and in what follows — is an integer. For this case 0** C Q for all 
E 


k. 
According to this decomposition, a free boundary I'*(r) is represented in 
(2 as the following set: 


nee aie Ts Ure) anes 
TP (r) = Ol ios T(r (r) = {aw € FF: |c(=)| = r(a, 1}. 


The structures of the pore space and the solid skeleton are characterized by 
elementary cells Y¢(r) Y,(r): 


Y;(r) = Int{y € Y : x(r,y) = 1}, Y,3(r) = Int{y € Y : x(r, y) = 0}, 


where Y = (—5) a) cR 

To the movement with increasing time of each component I*(r) of the 
boundary [*(r) in Q corresponds the movement of the boundary y(r) in Y: 
vr) ={ty EY: lyl =r}. 

As d,(r) in variables (y,t) we denote the velocity of the boundary 7(r) 
in Y in the direction of the normal n to y(r) in the pount y € y(r), outward 
to Yn): 

a(n) = -S"(@,1), (2.3) 
Similarly, as Dx,(r) we denote the velocity of the free boundary I(r) in the 
direction of the normal N* to I*(r) at the point « € I*(r), outward to 

Finally, for functions V of the type V = W(a,t, y) we define average by 
period 


(Wy = f vay, (L)y0) = [ xrwvay, (Vy) = [ (0-xtr.9)) vay. 
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In particular,for given structure r(a,t) with characteristic function y(r, y), 
given by formula (2.1), the function 


y) 


Wl rR 


4 
m(r) = i: x(r, y)dy = (x)y =1—- gar’ > 
Y 
is the Porosity of the solid skeleton in the point (x,t). 


2.1.4 Differential operators and matrices. 


Operator V without indices means differentiation with respect to the variable 
x; 
Operator V, means differentiation with respect to the variable y. 


1 
D(x,u) = 5 (Vou + (V,u)*). 
is a unit matrix. 


B:C=tr(B-C*), where B, C are Second-Order Tensors (matrices). 
a ® bis Dyad, for any vectors a, b, c: (a®b)-c=a(b-c). 


ji = 5(e @e! +e! @e’), {e!, e”, e*} is the standard Cartesian basis in 


R3 


A ®B is Forth-Order Tensor. 
(A @ B) : C = A(B: C) for any second order tensors A, B, C. 


2.2 Two — scale convergence [14]. 


In the present section we consider 1 — periodic in the variable y € Y = (0,1)? 
functions W(a,t, y) with (x,t) € Or. 


Definition 1. The sequence {w*} C La(Qr), is said to be two — scale con- 
vergent to the function W (x,t, y) € La(Qr x Y), which is 1 — periodic in the 


variable y € Y ( notation w* = W(a,t,y)), if for any smooth function 
ao =o(a,t,y), 1 — periodic in the variable y is valid the equality 


in ff ws(@.dotet =)de dt = 
[ Cf wet notet.adyaeae (2.4) 


The existence and basic properties of two — scale convergent sequences 
are proved in the following theorem: 


Theorem 1. (Nguetseng’s Teorem) 

1. Any bounded in Lo(Qr) sequence {w*} contains some subsequence 
two — scale convergent to some function 
W(a,t,y), W € La(Qr x Y), 1 — periodic in the variable y. 
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2. Let sequences {w*} and {eD(x, w*)} are uniformly bounded in La(Qr). 

Then there exists the function W = W(a,t,y), 1 — periodic in y, and 
the sequence {w*} such that W, V,W © Lo(Qr x Y), and sequences {w*} 
and {cD(x,w*)} (for simplicity we keep the same indices for subsequences) 
two — scale converge in L2(Qr) to W and D(y, W) correspondingly. 

3. Let sequences {w*} and {D(x,w*)} are bounded in L2(Qr). 

Then there exist functions w(a,t), w € Wy°(Qr) and W(«,t,y), W € 
Li(Qr x Y) A Wy(Y), subsequence from {D(x,w*)} such that the function 
W is 1 — periodic in y, D(x, w) € Le(Qr), Diy, W) € Lo(Qr x Y), and the 
sequence {D(x, w*)} two — scale converges to the function D(x, w)+D(y, W). 


Consequence 2.1. (Lemma B 13, Appendix B, [15}) 


ale 
Tf a(e)||[Vw*|loa, <M, where M does not depend on € and tim “© = 005 
E> 
then W(a,t, y) = w(a,t). 
Note, that weak and two — scale convergence are connected by the rela- 
tion: 


w*  W(a,t, y) > w*(a@) 3 | W(a,y)dy (weakly converges). 
Y 


2.3. Solenoidal functions in pore space. 
Let Y¢ (r(a,t)) =Int () Y;(r(a,t)). 
(#,t)EQr 
° 1 
By construction Y ¢4 @ because 0 # Y7(5) C Y;(r) for all 
1 
r(a,t), 0< r(x, t) < 9 


Recall, that some smooth function p(y) is called Solenoidal in vA y Evin 
We need the following statement, proved in [15] (Appendix B, Lemma B. 
15): 
Lemma 2.1. For any unit vector e there exists a solenoidal function p(y) 


° 1,0 ° O° 
such that p EW, (Ys), supp p CY ¢C Ye(r) and 


[ ewan =e. (2.5) 


2.4 Strong convergence criteria in L2((Qr). 


Definition 2. We say that the function u(a, t), bounded in L2(Q7), possesses 


a € Lo (0, T; Wz'(Q)), if 


ff pea mal ff |Vé|?dadt|2 
Qr Ot Qr 
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the time derivative 


for all functions € € W.''(Qr) with some positive constant M,, independent 
of €. 


Lemma 2.2. (Lions [26/, Aubin [27]) Let sequences {u*} and eae are 
uniformly bounded in the space L(Qr),and the sequence of derivatives erry } 
are uniformly bounded in the space L2(0,T; Wz" (Q)). 

Then there exists some subsequence of the sequence {u*} strongly conver- 
gent in Le(Qr). 

The generalization of this lemma for periodic structure with characteristic 

x 

function y*(x) = x() has been proved by Meirmanov, Zimin [28]: 


Lemma 2.3. Let y*(x) = vq \s where x(y) is 1 — periodic in y function, 
the sequences sa and {ve} are uniformly bounded in Lo(Qr), and the 
sequence {x* ~} is uniformly bounded in L2(0,T; W3'(Q)). 

Then there exists some subsequence of {c’} that converges strongly in 
L(Qr). 
Remark 1. We denote the norm of an element y in L2(0,T;W z'(Q)) as 
llc 


In this section, we prove a similar result for periodic structures with a 
special pore space structure. 


Theorem 2. Let the structure y(r,y) of the pore space be given by (2.1), 
where r € Mp and 


— a 
Mr ={r eH F Ap), 0 < r(#,t) < -, -6 < ar 


0<y<1, 0=const>0; Irlg, OVO May. 


(ety <0, 


Then any sequence {c°}, such that 


Iellaoy + IVEllany + Ie < M, 


where M does not depend on €, contains a subsequence, strongly convergent 
wn L2(Qr). 


Note that the non-positiveness of the time derivative of the function r € 
Mr means non-positiveness of the normal velocity of the boundary I*(r) 
in the direction of the normal N to this boundary, outward to the domain 

We divide the proof of the theorem into several steps. 
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Lemma 2.4. In conditions of Theorem 2 for almost all tg € (0,T) the se- 
quence {x °(a, to) C* (a, to)} converges weakly in L2(Q) to the function m(a, to) c(x, to). 


Proof. According to Theorem 1 the sequence {c*} (up to some subsequence) 
two — scale converges to some function c(a, t) in Le(Q-r): 


tim ff e*( me ake, =)dedt = Te yp) y dxdt (2.6) 
e0 Or Or 


for any test functions y € Le(Qr x Y). 

By definition of the time derivative of the function c* from the space 
L2((0, 7); Wz'(Q)) there exists some sequence of vector — functions {u*} 
such that ||w°||3.9,, < M* and 


a (at) oles) ae +u°- Vy) dxdt =0 (2.7) 
Qr 


for any test functions y, vanishing at the lateral boundary of the domain 
Or € R4, 

Here and throughout the text, MM denotes the constants depending only 
on Mo and T. 

Next we put y(a, t) = n(t) v(x). Then 


t= ff emermvaeat =i | [ e@.nnvacat (2.8) 


KO) = [ xe. vde, f(t) = f m(r)evde 


The limiting relation (2.8) means that 


Let 


T 


I=lim | a(t) fe (t)dt = / n(t) fo(t)dt. (2.9) 


e-0 0 


Coming back to the equality (2.7) with test functions of the form € = 
n(t) (a) we get 


i i (fs + nUS)at =0, Ut) = | uf (w,t)- Vu(w)de 


Q 


vi 
[ wrwrar< eo llVUllg < M* Valo. 


t) =UF(t), fe € W4(0,7), 


a) 
fS()|< My, [fol) — £5 (t2)| < My lt — ta]? 
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with some constant M,,, independent of e. 

The Arzela — Ascoli Theorem ([20] (Theorem 4, § 11, chapter II, p. 101) 
allows us to choose some subsequence { dots strongly convergent in the space 
C(Q) to some function f,(¢). 

Passing to the limit as <, — 0 in (2.9) {f;*} we arrive at the idewntity 


if T 
tim, f(t) r(tyat = f(t) Felt 
Ep,—0 ) 0 
for any test functions 7(t). 

That is, f,(t) coincides with f,(t) almost everywhere in (0,7), which 
proves Lemma. 


Lemma 2.5. In conditions of Theorem 2, there exists a subsequence {e;,} 
such that 


lim a | |V c** (a, to)|’dx = 0 (2.10) 
Q 


Ep—0 


for almost all to € (0,T). 


Proof. In fact, the uniform boundedness of the integrals 


II |Vc* (a, t)|°dxdt 
Or 


with respect to ¢, implies equality 


lim af | |V c** (a, t)|/°dadt = 0. (2.11) 
Qr 


ER,-0 


Let 
Ux(to) = é | | Veo" (a, to) Pda. 
Q 


Then (2.11) means, that the sequence {u,} converges strongly to zero in the 
space L,(0, 7’). 

Referring to Theorem 1 ({20], Chapter 37, § 1, p. 379), we conclude that 
there is a subsequence {ux}, which converges to zero almost everywhere in 
(0,7) as k > oo. 


Lemma 2.6. In conditions of Theorem 2 any bounded in the space ILo(Q) 
sequence {C**(a,to)} converges weakly and two — scale in the space Lg(Q) 
(0,7) to some function c(a, to). 


Proof. Let the sequence ¢; be the same, as in Lemma 2.5 and II C (0,7) be 
the set of full measure for which the condition (2.10) holds true. 


16 


The boundedness of the sequence {C**(a,to)} in the space Lo(Q) im- 
plies the existence of some sequence (we leave for simplicity the same in- 
dices), which two — scale converges to 1 — periodic in the variable y function 
C(x, to, y) from the space L2(Q x Y). 

Integration by parts of the expression ¢,V C**(a, to) - v(~)u(a) results 


Ek 
the identity 
zs z 
cx f VE*(a to) @(=)o(@)de = 
Q Ek 


— 24 f (0, to)(ol 


x a x 

2). Vu(e))de— f e(«,t0)(V (=) Hle)ae, 
Ek Q Ek 

which holds true for arbitrary functions yp € WAY) and w EW (Q). 


Passing to the limit as ¢; — 0 we get 


i, (a) ( . Cle, to, ¥)V «ely)dy) de = 0. 


Due to the arbitrary choice of functions y and w, the last identity is equiv- 
alent to the equality C(x, to, y) = c(a, to). 


Lemma 2.7. In conditions of Theorem 2, the sequence {c**(x,t)} strongly 
converges in Lo(Qr) to the function c(a,t). 


Proof. To prove Lemma we put 


H' =W", (Q) C H® =L,(Q) CH! =W~*, (Q), wy (aw, to) = C* (aw, to)—c(@, to). 


Taking into account inequality (9) ([29], § 12, Chapter I, p. 50) we obtain 


l|we(-, to) [liso < 7 [lve (-, to) lian + » lee, to) Ilr. 


Next we integrate the last inequality with respect to time 


Si 


T T 
i lue(.2)eat <n f lun. HIKadt +e, f l|we(-, t)lieadt. (2.12) 
0 0 0 


and after that use a compact embedding of the space H°(Q) into the space 
H~'(Q) ({30], Theorem 3, § 2, Chapter 4): the weak convergence of the se- 
quence {w,(.,t)} in the space H°(Q) implies the strong convergence this se- 
quence in the space H~'(Q). 


That is, jim \|wz(., t)|[f-1dt = 0. 
00 
The last equality (2.12) and arbitrary choice of the positive number 7 


- 
means that iim f \|we(., t)||zodt = 0. 
k-0o 0 
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2.5 Extension Lemmas. 


The results on the extension of solutions of differential equations from one 
domain to another one are very important in the theory of homogenization 
of differential equations (see Zhikov et al. [10], Zhikov{11], Conca[31], Acerbi 
et al. [32]). 

For example, if a sequence has different differential properties in different 
domains, and when homogenize it is desirable to preserve the best differential 
properties of solutions, then the better choice would be to extend the solution 
from the domain where it has the best differential properties onto the domain 
selected for homogenization. 

All the results known up to now have been obtained for periodic struc- 
tures. Fortunately, 

the results of Lemma 1 ({10], §1 of Chapter III, p. 95) on ”soft inclusions” 
are fully applicable to the case of the non-periodic structure with special 
periodicity of the pore space under consideration, when this structure is 
determined by the relations (2.1). 

The structure that we consider here is, in principal, the ” soft inclusion” 
and the special periodicity of the pore space does not play any role here 
because the extension is done for each individual cell Q** © across the 
connected component of the boundary P'**(r) = T*(r)NQ** of the boundary 
Pe{¥). 

The following statement is true 


Lemma 2.8. Let r € Np and functions r(a,t) define the structure x(r,y) 
of the domain QF 7p, given by relations (2.1). 


Then for any bounded in the space W3”°(Q5.r(r)) sequence {c°} there 
exists an extension from the domain Q;7(r) onto domain Qyr such that 


& € W)°(Or), [le llaar <M lle*lla.0¢ 0 
I[Ve‘lloa7 < M ||Ve*lla.a¢.(7): (2-13) 
Next, we will need functions that take the given values on the connected 
components I'**(r) of the boundary ['*(r). Namely, we will need an exten- 


sion of the vector function v* (a,t) from the boundary [*(r), defined by 
function r(a,t) by equalities (2.1), onto domain Qr. 
At the same time, the connected components of the support of this func- 
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tion will be subsets of set 0**: 


ve o ke o ke o 1,0 
O° (eat) vy (r;a,t), v EW, (OF”); 
keZ 


IV lo, +¢||D(z,v* ) lla, < ¢M Mo; 
V- v¢ (r; x,t) =0, (@,t) € OFF, 

8 x 5 

*(r32,t) =0 —)|>—; 

Vv (r; a, ) Is) 2 753 

v* (r;X0,t) — ( vw (Cima t) N*)N® =f), 

v® (r3a,t)-N® = —N* D£((ao,t)), ao ET **(r), (2.14) 
where N* is the unit normal to the boundary ['*(r) at the point a of this 
boundary at the time moment ¢, outward with respect to the domain 0;(r) 


and Dx, is the velocity of the boundary [’*(r) in the direction of the vector 
N°. 
Lemma 2.9. Let r € Nr and functions r(a,t) define the structure x(r, y) 
of the domain QF p, given by relations (2.1). 

Then for alle > 0 there exist functions v* (r; x,t), satisfying the condi- 
tions (2.14). 

1 
Proof. Let r° = eae rw,t) < 5 and infinitely smooth periodic function 
xz tyEeQp 


So(y) be selected from the conditions 


acta ra lul <r, 


0 ly| > 3 


We consider a connected component ['°*(r) of the boundary ['*(r), where 
conditions (2.14): 


x 1° r 
Iso(=) = r(o,t), 5 v* (7; Xo, t) = Diy (#0, t) so(=), ay €T**(r) 


(2.15) 
must be fulfilled. 


We look for an extension v* (r;a,t) of the function v* (r; a, t) as 


° 


x x 
ve (r; a, ¢) = —6 Diy (@, t) So(=) ue(@, 2) 
where &p = Ek+e (2), w%=ek+e s() and the function to be find should 


x 
be u;(a, = 
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We have to extend this function v* (r;ao,t) from the boundary ['**(r) 
into a ” liquid component” of the domain Q** so that the extended function 


ve (r;a,t) should be solenoidal. 


The required solenoidality of the function v* (r;a,t) provides the partial 
differential equation of the first order 


02>" (at) = 


E 
: x x 1 x x 
Dy (2, 1)( (so(=) Vn Ur(@, =) oe (so(=) Vy u(x, =))) (2.16) 
for the function u;(x, a completed with initial condition 
u(@o,—) = 1 (2.17) 
at the surface 
Per) = {ag EO" 55 =ek+s(—)}. (2.18) 


For a while, we will omit the upper index e. 
We rewrite equation (2.16) as 


a, (2, =) Vat + = ay(2. =) Vy uz = be(ax, =) ur(a, 2 (2.19) 
where 
ba, 2) = —( (Vz Dw(w,t) -60(=)) + = Dw(w,t) Vy -s0(=)); 
a,(x, =) = Dy(z,t) so(=), (x,t) € Ap. (2.20) 


x 
Our immediate goal is to extend the function u;(ao, =) from the surface 


[**(r), given by the relations (2.18), into the liquid part ee) of the 
domain ** along characteristics of the partial differential equation of the 
first order (2.19). 

Let s be the extension parameter, 2 be the starting point on the surface 
[**(r) for the characteristic of equation (2.19). 

According to [33], the Cauchy problem (2.17), (2.19) is equivalent to the 
Cauchy problem for the system of ordinary differential equations 


Ox 1 
De 6 Xo) = at (x(s, Xo), 22s, Zo)), x(0, Lo) = Lo; 
1 
Fe (8 Xo) = by (a(s, Xo), e x(5, ge 
U,(0, #9) = 1 for Is(—) =r(ao,t) (2.21) 
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1 
for new unknown functions x(s, a) and U;(s, xo) = u(a(s, xo), Z a(S, @o)). 
Due to the uniform boundedness of the right hand — sides a;(a,— a) 
E 


1 Za qe 
and b.(a, x) in the space C!(Y x Or), the Cauchy problem (2.21) has a 
unique solution {x(s, a), U:(s,@o)} for all (a,t) € Qr such that U;(s, ao) 


and —(s,ap) are continuous functions on any bounded interval s € [0, S] 


[34] where 
1 
Iso(= a(S, ao))| > 0. (2.22) 
A natural question arises where the point x(s, a9) will start to move from 
the initial position with increasing parameter s? Into a liquid or a solid 
component? Since the vector ¢9(—) at the boundary point directed into the 


liquid component, then for s > 0 the point a(s, ag) will begin to move into 
the liquid component with the growth of the parameter s: x(s,a) € OFF (t) 
for s > 0. 


1 
Direct differentiation of equality U;(s, ao) = uz(x(s, @o), = z(s,@)) with 


respect the parameter s gives us one more equality 


OU, 1 6) 
—*(s, a9) = Va (a(s, ao), - a(s, ao) - ele Lo)+ 


Os Os 


1 1 Ox 
is Vy uz (x(s, &o), Z a(s, a) : De (8 #0) = 


Vz ur(@(s, 20), . x(s,&o)) - a:(x(s, Xo), ~2(s, &o))+ 


1 1 Ox 1 
. V, u(x(s, £0), A x(s,@o)) - race Zo) + ax(x(s, Lo), =a, @o)). (2.23) 


Thus, the second equation in (2.21) and equation (2.23) prove the validity of 


the relation (2.16)and, that is, the solenoidality of the function v* (r; a, t) 
in the domain OF) everywhere wherever condition (2.22) is satisfied. 


On the other hand-side, where this condition fails, one has that v* 
(r;a,t) = 0. So, indicated function will be solenoidal throughout the liq- 
uid component. 


Since throughout the text of the article, the values of the function v 2 
(r;a,t) are used only in the liquid component, the extension of this function 
into the solid component is standard. For example, along the normal to the 
boundary with further multiplication by a cutoff function. 
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2.6 The Leray - Schauder fixed point theorem. 


Definition 3. A continuous operator T acting from the Banach space X into 
itself (notation: T : X + X) is called Completely continuous if it maps 
every closed bounded set into a compact one. 


Theorem 3. Let the completely continuous operator T(X, x), acting from the 
bounded Banach space X of elements x € X into itself, continuously depends 
on the real parameter \ € [0,1] and let for \ = 0 the operator T(0,x) there 
is a fixed point x9: T(0, 20) = 20. 

Then for all X € [0,1], the operator T has at least one fixed point x), such 
that T(A, xy) = 2). 


3 The problem statement. 

Let the function r*(a,t) determine the structure of liquid and solid compo- 
nents in the domain Qr in according to the formulas (2.2), where ro(a) be 
1 Pe 
a given function, 0 < r*(a,t) < 5 and r*(ax,t) = max{0, ro(a) — r*(a, t)}, 


and 


, p= — (3.1) 


are dimensionless variables. 
We look for the solution {7*, v*, p°, c°, G}, j = 1,2,...,n to the system 
of differential equations 


V - (af D(z, v*) — pT) = 0, V-v*® =0, (3.2) 
(elon 
an V-(doVeE —v* &), (3.3) 
Oe 
ai, +V- (vc) =0, j=1,...,4, (3.4) 


in unknown domain 05.7(7°) describing the dynamics of the interaction of 
an acid solution in a carrier fluid with a concentration c*(a,t) with a solid 
skeleton 05 (7°). 

In (3.2) — (3.4) pg — atmosphere pressure, v* and p* — velocity and pressure 
in the liquid, c; — concentrations of products of chemical reactions at the free 
boundary, P'*(7*) = 005(7*) NOQZ(r*). 

Generally speaking, the dimensionless diffusion coefficient dy can depend 
on the concentration of the acid. In our setting, we consider it to be a given 
positive constant. 

At the free boundary [*(7*) hold true boundary condition (1.1) and con- 
ditions 

Un = —6 Dy, (3.5) 
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ve =i N*, (3.6) 


Oc 
(Dy — Un) & + do aN =r e; (3.7) 
Co eno ees (3.8) 


expressing the laws of conservation of mass and momentum [15] — [17]. 


Here viy ay. N® is a normal component of the liquid velocity v*, 
Oc® 
aN {Ve-.N ©, N® is a vector of the unit normal to the free boundary 


P*(r*), outward to the domain Q;(7r*), 6 = eee Ps and py are densities 
of solid and liquid components respectively, 6° and 65, j = 1,...,n are given 
positive constants. 

Finally, at the given boundaries with injection wells S' and producing 
wells S?, and at the impenetrable boundary S°, the following conditions 


(af, D(x, v*(a, t)) — p°(a, t) ) n(x) = 


—pn(x), ©€ $7, t>0, 7 =1,2, (3.9) 
Oce 
=0, 2:6 8°, £>0, 3.10 
an a (3.10) 
v =0, x€ S°, t>0, (3.11) 
(x,t) =c°(x), «Ee S'US?, t>0, (3.12) 
Clea) 0) Tce Gre Se SO (3.13) 


are met. 

In (3.7) — (3.13) n is the unit normal vector to the boundaries $°, S' and 
cigs 

The problem (1.1), (3.1) — (3.13) is completed with initial conditions 


Ce. =en) cae =0, 7 = lan, 2 eQ, (3.14) 


F*(w,0)=0, EQ. (3.15) 


Note, that the problem (1.1), (3.2), (3.3), (3.5) — (3.7), (3.9) — (3.12), (3.14), 
(3.15) on determining functions {7*, v°, p*, c} (let call it as Problem A*) is 
solved independently of the problem on determining functions cj(x,t), 7 = 
12, ces Fb: 

The latter are determined after finding the functions {v*, p*, c*}. 

1 
We will assume that p' = const, i = 1,2 and p°(x1) = p'(a1+ 5) +p" (21 — 
). 


1 
Z 
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For new unknown functions v°(a,¢) and p*(a,t) = p*(a, t) — p°(1) equa- 
tions (3.2) in the domain 0p will take the form 


V- (of D(a, v°) — pl) = Vp, (x,t) € 2% 7(F*), (3.16) 
V-v° =0, (#,t) € OF r(r°), (3:17) 

and boundary condition (3.9) will be rewritten as 
(oo Diao) pl) oN =0,4 € 83,7 = 1,2; > 0: (3.18) 


Let us transform the equation (3.3) and the boundary condition (3.7) so that 
the transformed equation and the boundary condition are equivalent to some 
integral identity. 

Namely, we put 


Oc® Se E & Be € rE 
a =" (doVci — v° (c Pee) (x,t) € OF 7(F*), (3.19) 
Oc® E E Be = cf e 
do ON —UyN (c + re = 0, (a6) E TB(r a (3.20) 


By its original definition, the acid concentration c* is always nonnegative, but 
its nonnegativity in our model does not follow from anywhere, since for the 
equation (3.19) and the boundary condition (3.20) there is no the maximum 
principle. 

Therefore, we temporarily change the equation (3.19) to equation 


Oc® 


ae V- (doVet — v' W(c%)), # € N4(F*), t > 0, (3.21) 
and the boundary condition (3.20) to the boundary condition 
Oc® a 
dos — yi ule) =—c® Dy, # ET(F*), t> 0. (3.22) 


In (3.22) N is a unit vector to the boundary I*(7*) of outward to the domain 
O-(7*) and 

v(s)=94+ 5, 0s <0 vis) =). 620; sac: (3:23) 
We call the equation (3.21) as Modified diffusion — convection equation 
for the acid concentration, and the boundary condition (3.22) as Modified 
boundary condition for the acid concentration. 

By problem Problem B*(r) we mean problem A*® without boundary con- 
dition (1.1), but in a given domain 2;(r) with a pore space structure x(r, y), 
is determined by the function r(a,t) in accordance with the relations (2.1). 

Finally, by Dynamic problem B*(r) we mean initial — boundary value 
problem (3.11), (3.16) — (3.18), by Diffusion — convection problem 
B*(r) we mean initial — boundary value problem (3.10), (3.12), (3.14), (3.19), 
(3.20), and by Modified diffusion — convection problem B*(r) we mean 
initial — boundary value problem (3.10), (3.12), (3.14), (3.21), (3.22). 
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4 Equivalent integral identities. 
Formal homogenization. 


To formulate main result we need to know homogenized models H and H(r). 
To get the differential equations of macroscopic mathematical models, we 
must somehow homogenize the microscopic mathematical models A* and 
BS(r). 

Due to the remarks made earlier, it will be enough to homogenize only 
the model B*(r). 

We assume that 


FEEM, ve, ch © We°(DEr(r*)) FF ELo(Q8y), |p°|g? = Mo < 0. 


Following the known schema [15] we extend solutions v* and c* onto the 
domain 9% -(r°) and then write down the problem A* in the equivalent form 
of the corresponding integral identities so that the conditions on the free 
boundary will be included in the integral identities. 

According to the lemma 2.9, the function c® is an extension of a function 
c®, and the Lemma 2.8 was devoted to the extension of the function v€. 

By construction (see Lemma 2.9) the function v- belongs to the space 
W,°(Qr) and satisfies to the boundary conditions (3.4). 

Now we put 


5*(@,t) = eae as (#,1) € %r(F*), 
; v (F&:a@,t) as (x,t) ENE ,(F*) 
= _ |p (@,t) as (w,t) € 0F7(7*), 
ee i as (x,t) € 0 oF). ey 


Then the function 0° € Wy°(Qr) is an extension of the function v* from the 
domain 04 7(7*) onto the domain Q{7(7*) and satisfies boundary condition 
(At). 


Definition 4. Let the structures y°(a,t) = x(r(a, t), =) of the pore space 
OQ; 7(r) be defined by a function r € Mr. 

Then the solenoidal function 0° € W.°(Qr) and functions p* € L2(Qr) 
and ¢° € Wy°(Q,) are called the Weak solution to the problem A‘, if the 
integral identities 


al x* (a%, D(x, ¥°) — p*1) : D(z, y) dxdt = 
Or 


-| ve -paeat— ff ax V% F :D(x,)) det, (4.2) 
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es (3% 45°. Wh) drat =0. (4.3) 


Fe (x, 0) c* (a, 0) €(a, 0)da+ 


i [x(- a + VE. (iy Ve — 8 (@ +) andr = 0 (4.4) 


and boundary conditions (1.1) and (3.11) are fulfilled. 

The identity (4.2) holds for all smooth functions y, vanishing at the 
boundary $° of the domain 2 for t > 0. 

The identity (4.3) holds for all smooth functions ~, vanishing at the 
boundary S' U S? of the domain Q for t > 0 and at the initial moment of 
time. 

The identity (4.4) holds for all smooth functions €, vanishing at the 
boundary 5! U S? of the domain Q for t > 0. 


Remark 2. A little below we will show that the boundary condition (1.1) is 
equivalent to the is equivalent to the integral identity (4.25). 


The identity (4.2) is equivalent to the dynamics Stokes equation (3.16). 
In this identity 


Oo € ° oO € Oo — 


P= /aiD(e,v°) F €Lo(Qr) || F 


laa, < MM, 


where the function v* is defined by the formula (2.14). 

The identity (4.3) is equivalent to the solenoidality of the function v * (a, t) 
and boundary conditions (3.5) and (3.6). 

The identity (4.4) is equivalent to the diffusion — convection equation 
(3.21) and boundary and initial conditions (3.7), (3.10), (3.12), (3.14) and 
(3.22). 

The equivalence of differential equations with corresponding boundary 
conditions to integral identities follows from the formula of integration by 
parts ({16], part I, §12) in the form 


Pee e(Sot+v- Bydoar + f° as S 4B. VE)dxdt = 
2; (6) Qf (t) 


i E(x, to) A(x, to)dx -f &(a, 0) A(x, 0)da+ 
5 (t) 


0 f(0) 


dub &(B.N-—ADw) singdodt, (4.5) 
0 


Here N € R? is the unit normal vector to I(t), pointing outward to 0,(t), 
Dy is the normal velocity of the boundary I(t) in the direction of the normal 
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N, and yg is the angle between unit vector J of the time axis and unit normal 
vector v € R* tor, pointing outward to the domain 2,7, such that sin y = 
v-Nandcosy=v-l. 
For example, in the identity (4.4) A = ¢* and B = v* (c+ “do Wie: 
a 
The transformation of the diffusion — convection equation (3.3) leads to 
the equalities 


E 


aed 


to ag 7 Be 
—-—&+VE-(do Vee —v* (CF + dxdt = 
i a ot : ( : ( a), 
to ie : Bae Mews ie OGE 
: [8 Dy + Un (CF + ee — do aN? sin y dodt+ 
to O€ é ae ore GE 
| bleu + VE- (do VE —0* (E + —))) deat. 


2 Sopcast 1G Oe. ? 

Integrand expression Vg = —C* Dy +0xy (C© +=) —do ce in the integral over 
the free boundary is equal to zero due to the boundary conditions at the free 
boundary: 

€ EOE. an iE ETE Oc® ero. eet oi eG. 

Dy =a&c®, Un = —6 Dy = —d0° C*, — do ae =F c+ DNC — UNC. 

Definition 5. Let the structure y(r,y) of the pore space 0; 7(r) be given 
by the function r € Np. 

Then functions v* and p* are called the weak solution to the Dynamic 
problem B*(r) if integral identities (4.2) and (4.3) and the boundary condi- 
tion (3.9)are fulfilled. 


Definition 6. Let the structure y(r, y) of the pore space 0; 7(r) be given by 
the function r € Np and function v~ be the weak solution to the dynamic 
problem B*(r). 

Then the function c* € V2(Q7) is called the weak solution to the Diffu- 
sion — convection problem B*(r) if the integral identity (4.4) and bound- 
ary conditions (3.10) and (3.12) are fulfilled. 


Definition 7. Let the structure y(r, y) of the pore space 07 7(r) be given by 
the function r € tp and function v* be the weak solution to the dynamic 
problem B*(r). 


vf 


Then the function c* € Vo(Q,) is called the weak solution to the Mod- 
ified diffusion — convection problem B*(r) if the integral identity 


[ ee - oS + ve. (do Ve" — B°w(@*)) )deat+ 
[xi eto)e (ato) gle, to)de— f (2, 0)e"(e0):é(e, 0)da = 0 (46) 


Q 


and boundary conditions (3.10)and (3.12) are fulfilled for all to, 0 < to < T 
and for all smooth functions €, vanishing at the boundary S°. 


(4.6) u(s) (3.23). 


Definition 8. Let the structure y(7, y) of the pore space 0; 7(r) be given by 
the function r € Ip and function v~ be the weak solution to the dynamic 
problem B*(r),and the function c* be the weak solution to the diffusion — 
convection problem B*(r). 

Then the triple of functions {v°, p*, c*} is called the Weak solution 
to the problem B*(r). 


Lemma 4.1. Let the structure x(r,y) of the pore space Q7_(r) be given 
by the function r € Mr and the solution {v*, p*, c*} to the problem B*(r) 
exists and has the smoothness necessary for the selection of convergent sub- 
sequences: 


p, v°, eD(a, 0°) €L2(Qz), & € Wy (Qr), 


ae € L2 (0,7; Wz*(Q)); 
IP“ llo0r + IP *ll20% + ILVou D(z, 0°) Il2,07-+ 
po Oc * 
IV Fler tls lw <M Mo. (47) 


As M here here and everywhere else, we will denote constants that do not 
depend on €. 

Then there exist some subsequences of the sequences {v°}, {p*} and {c*} 
such that 


B°(a,t) > v(a,t), B°(a,t) "> V(a,t,y), 
eD(x, 0° (a, t)) = D(y, V(a,t, y)); 
D(a, t) — plat), P°(a,t) "> Py (a, ty) x(a, ty); 
ve (x,t) strongly converges in Lo(Qr) to zero; 


oe 
F (x,t) strongly converges in Lo(Qr) to zero; 


C*(a,t) strongly converges in La(Qr) to the function c(a, t); 


D(x, é*) “> D(z, c) + D(y, C(z,t, y)) 
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as é—+ 0. 
Here 
v €LA(Q7), c, pe Wy°(Qr), Ty, V, C, Diy, V), D(y, C) € Le(Qr x Y). 


The proof of the Lemma follows from the results of § 2 on the choice of 
convergent subsequences. For the limiting procedure, it is necessary to know 
the behavior of the coefficients a/,, 

The limit as ¢ > 0 in (4.3) is quite standard (see [15], pp. 6 — 11): 


1 
v=——C(r) V(p+p"), (@,t) € Or. oe 
1 
om 
VSO (a, t) € Or. (4.9) 


At the same time identities (4.2) and (4.3) contain boundary conditions 
in =v-n=0, ES, p=0, eES'US?, t>0. (4.10) 


on the given boundaries. 
The limit in the identity (4.4) as « — 0 gives us the identity 


[ EP Va ee de [ HE DOC) O aes 


Lf cmeg$ + ve: oe) Ver wle+ 9) )dedr = 0, (4.11) 


which holds true for all smooth functions €, vanishing on the boundary S'US? 
of the domain (22 for t > 0. 
The identity (4.17) is formally equivalent to the Homogenized diffu- 
sion — convection equation 
O ‘ 60 
ne) = WV: (do C*(r) - Ve v (c+ 5). (4.12) 
BF and af ase — 0. 
The constants G° will be defined in (4.23). 
In the present publication we restrict ourself with the case 


E 


[i =] 0 0 a << oo, ine =0'= ‘const (4.13) 


0 aé 


of Weakly viscous liquid (sce {15]). 
The limit as ¢ — 0 in (4.3) is quite standard (see [15], pp. 6 — 11): 


o= =r) -V(ptp’), (x,t) € Qr. (4.14) 
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Viv a6, (x, t) € Qe. (4.15) 


At the same time identities (4.2) and (4.3) contain boundary conditions 
Un =v-n=0, ES, p=0, xe S'US", t>0. (4.16) 


on the given boundaries. 
The limit in the identity (4.4) as « — 0 gives us the identity 


[ mea.to en ce? (a) E(x, 0)dx+ 


ie yA — mes + VE. (do C°(r) -V Ve-v(e+5)))drdr = 0, (4.17) 


which holds true for all smooth functions €, vanishing on the boundary S'US? 
of the domain 22 for t > 0. 

The identity (4.17) is formally equivalent to the Homogenized diffu- 
sion — convection equation 


a ; i 
Hvme) = V: (doC(r)- Ve—v(e+ 5)). (4.18) 


Symmetric strictly positive definite matrices C’(r) and C°(r) depend on the 
structure of the pore space Y;(r) and are defined by formulas (1.1.27) and 
(10.1.61) in [15] accordingly: 


can f vie (r; y) ) @e')dy, 
Yp(r 


re —V,I' = -e', V,-V'=0, |y|>r, 
Vitrsy) =0, lyl=r, (4-19) 


(Vie t+e')-n=0, |y|=r, (4.20) 


where {e', e?, e®} is standard orthonormal basis in R?. 
Diferential equations (4.14), (4.15) and (4.18) are completed with bound- 
ary conditions (4.16) and boundary and initial conditions 
Oc 0 


Fy WO = 0, BES, cHe, BESTUS’, t>0, (4.21) 
n 


30 


c(w,0)=c°, EN. (4.22) 


We call as the Problem H(r) the initial boundary — value problem (4.14) 
— (4.22) for the given structure y(x,t, y) of the pore space, defined by the 
function r € Mp. 

The boundary — value problem (4.14) — (4.16) we call as the Dynamic 
problem H(r). 

Finally, the initial boundary — value (4.18) (4.20) — (4.22) we call as 
Diffusion — convection problem H(r). 

Everywhere below, where it does not cause discrepancies, we will write 
C’ and C* instead of C’(r) and C“(r). 

When deriving the problem H(r) we did not use anywhere an additional 
boundary condition (1.1) on the free boundary I(r). 

That is, the problem H(r) is a homogenization of the problem B*(r). 

Thus, to derive the Problem H as a homogenization of the problem A* we 
simply complete the problem H(r) with condition, which is a homogenization 
of the additional boundary condition (1.1) on the free boundary. After that 
we find the structure x*° of the pore space according to the formulas(2.1). 

To homogenize boundary condition (1.1) we will need the following 


Assumption 1. For the problem A* hold true conditions 
ao =e), 6° Ss, (4.23) 
where @ is a given positive constant. 


Lemma 4.2. Let functions c*(a,t) be extensions of functions c*(a,t) (see. 
Lemma 2.8), the sequence {&*} be bounded in the space W3°(Q) and lim \|e*— 
E> 
Cllo.ar = 0. 
Then under condition of the Assumption 1 the equality 
d,(x,t) =@c(az,t), ye y(a,t) CY. (4.24) 
will be a homogenization of the additional boundary condition (1.1) 


Proof. The additional boundary condition (??) on the boundary [<(r) is 
equivalent to the integral identity 


i x8(—S(a®-€§) +eV- (C25) dat = 0 (4.25) 


x 

which is valid for any smooth functions €7(r,a) = &.(r, =), functions ¢, 

vanishing at t = 0 and at t = tg and at boundary OQ, and functions 
x 

a.(r,x) = a.(r, =), such that a, vanishes outside of some small neighbour- 


hood of y-(r) and a.(r,y) = n-(r), where n,(r) is the unit normal to 7,(r), 
outward to the domain Yf(r). 


dl 


In fact, let ¢(a,t) be an arbitrary function, vanishing at t = 0 and at 
t = to. Then, according to the identity (4.5) of integration by parts with 
A=1and B= 0, we have 


_ bow oe dt = -{ Dy u sing do dt. 
€ T(r) 


In the last identityy is the angle between the time axis and the unit normal 
N* to the boundary I'*(r) outward with respect to 25(r), Dy is the velocity 
of I'*(r) in the direction N°. 

Next in the integral over the surface [7(r) we use boundary condition 
(1.1) and Assumption 1: 


ae. andl = -{{ a CFU singdodt =~ ff eOc* u sin y dodt. 
29 to (r) to (r) 


Let now 
u=a‘(w,t)-€°(@,t), a°(w,t) = a(w,t, =), €°(@,t) = ¢() &o(=, 2). 


with 1 — periodic in y functions a(a,t, y) and €)(y,t) such that a(a,t, y) = 


n(x, t,y) for y € 7(a, t). 
One has 


yy 


an c (a°-&° save + [ J e0c*(n- €*) sinydodt = 
E(r) 


to 


ee <a" -€5)+e0V- (¢ & €5))dadt = 0. 
The obtained identity 
a x(Z(Gat- 65) +e0V- (C25) drat = 0 (4.26) 


is equivalent to the boundary condition (1.1). 
x 
It is valid for any functions a*(x,t) = a(a,t, Z) , a(x,t,y) = n(a,t,y), 
where n(a,t, y)is a unit normal vector to y(a,t), outward with respect to 
the domain Y,;(a,¢), arbitrary functions €)(y,t), equal to zero at t = 0 and 
at t = to, and arbitrary functions ¢(a), vanishing on the boundary OQ. 
Passing to the limit as « > 0 in (4.26) we get 


feo LL, (F(a &) + Vy (@e&s)) dydt = 0. 
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Here we have used the last statement of the Theorem 1 on a two — scale 
convergence of functions from W,°(Qr) to the function independent of the 
fast variable y. 

The equality (4.24) we need follows from the last identity for the arbitrary 
functions ¢ and from the identity (4.5): 


0= [ [ (5(@- 0) + ¥y (obs) dt - 
. 


is [ra +6c)(n-&,) sin ydodt = 0. 
0 vy 


Turning to the boundary condition (4.24) we see that its fulfillment entails 
equality 


The last relation suggests a scheme for proving the existence of a solution to 
the problem H using Schauder’s fixed point theorem. 

We fix the set 27 of functions r(x, t), defining the pore space (7) in 
the variables (a#,t) and the pore space Y;(r) in the variables (y, f). 

For given structure of the pore space we find the solution 

{v =F’(r), Vp = F?(r), c= F“(r)} to the problem H(r) and function 


R(x, to) = ro(a) — 6 / c(a,t)dt = F(r), (4.27) 


by which a new structure of the pore space is constructed in accordance with 
the formulas (2.2). 

The fixed points r,(a,t) of the operator F(r) define the characteristic 
function x(r., y) of the pore space Y;(r,) according to the formula (2.2) § 2, 
for which the problem H(r,) coincides with the problem H. 


Lemma 4.3. Let 


Om 
Vm, OE — ee Cee Be 


plz? < (C?-2)-2< pleP, wale? < (C-2)-2< ple? 


and p, © € Lee(Qr) NW5°(Qr), v € Lao(Qr) are the weak solution to the 
problem H(r). 
Then operators F°(r), F’(r) and F*(r) are defined correctly. 
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Proof. The difference p = p,; — po of two possible weak solutions p; and po 
of the problem (4.14) — (4.16) satisfies in the domain 2;, the homogeneous 
elliptic equation 

V-(C’- Vp) =0 (4.28) 
and homogeneous boundary conditions (4.16)at the boundary S of the do- 
main 2 for all t, 0 <t < to. 

Strict positive definiteness of the matrix C’, guarantees the uniqueness of 
the solution of the problem (4.16), (4.28) and correct definition of operators 
F’(r) and F?(r). 

Now consider the difference u = m,c, — M22, m; = m(r;i), 7 = 1,2, 
two possible solutions cy and cp to the problem (4.18), (4.21), (4.22) ,which 
satisfies the difference of integral identities (4.4): 


[ wet) €(a, to)dx+ 
Q 
[ [ews a ve CoV ijdedt = 
” Finalise 
a mom 


where € is arbitrary smooth function, vanishing at 2 € S'US? and 0 < t < to 
and the following homogeneous boundary and initial conditions 


u(z,t)=0, wnES'US?, O<t<T, 


O 
ay =O xe S°,0<t<T, u(x,0) =0, x EN. 


Using the standard time smoothing procedure ([{25], chapter I, § 4)with the 
test function € = (u,); and inequality ab < \a? + C) b? we get 


ge aa 


[[? (un)z,°(C°-Vu)dadt < | |(Vu) 7, Pande +Cy ff u?dadr. 
Q, ™M er Q¢ 


The limit as h > 0, strict positive definiteness of the matrix C°: 


e(-(C*-€)) > we? 


the choice is small enough \, give us a differential inequality 


[ews dx <M? II u’(x, t) dxdr.u(a, 0) = 0, 
Q on 


Gronwall ’s lemma (Lemma 5.5, § 5,chapter II, [25]) implies c(a, t) = 0 almost 
everywhere in 7 and, that is, the correct definition of the operator F“(r). 
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5 Main result. 
Theorem 4. Let c° € C?*7(Q), and the conditions (4.13) and condition 
0< (a) < ch =const <1 EN, t>0 (5.1) 


be fulfilled. 
Then for all r € Mp the problem B(r) has a unique weak solution 
{v*, p°, c©} such that the estimates 


IP llaog pir) + [lO llao¢ 0) + [Ven D(a, v*)la,.02 07) S 


M(Mo)|V pl? = Mp, (5.2) 
llellaag.@) +1V eFllaag 0) < M(Mo) (llS*? +.M,)=M., (5.3) 


Oc © 
Se ivy CNY, 5.4 
1 lw < MM.) (5.4 


O<c*(a,t) <c*, (x,t) EOF7(r) (5.5) 
hold true. 
Here constants M(M,) and M(M_) do not depend on «. 


Theorem 5. Let in conditions of Theorem 4 be fulfilled the Assumption (1). 
Then for all T > 0 the problem H(r) has a unique classical solution 


{v, p, c} such that Vp, v € HY? (Qr), ce He+1 3" (Op), 


[V pl + fol < MM) (5.6) 

and 
0 < c(a,t) < c*, (x,t) € Or, (5.7) 
Ilr? < M(M.). (5.8) 


Theorem 6. In conditions of Theoremd for all T > 0 the problem H_has a 

unique classical solution {r*, v, p, c}, such thatr* € Nr, Vp, v € H 13 (Or), 
24+7 = 

ce H2+1 3" (Gp) and hold true estimates (5.6) — (5.8). 


6 Proof of Theorem 4. 


According to [25] (§ 9, chapter IV — correctness of the diffusion — convection 
problem) and [35] (Theorem 2, § 5, chapter III — correctness of the dynamic 
problem), 


é& ents (Qer(r)), v, Vp’ Ela (0,7; W? (24(r))) NH? (Q;7(r)), 
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where g > 3(1+7) is any. 

Note that these theorems are formulated for cylindrical domains. But 
with the help of local estimates, all statements can be proved for non- 
cylindrical domains. 

The belonging of the function v* to the space H™? (Q;-r(r)) is proved in 
the same way as in the Lemma 7.5 (estimate (7.25)). 

The dynamic problem B*(r) is a linear one. Therefore, to prove the 
existence of weak solutions to the problem B*(r), it is sufficient to obtain 
appropriate a priori estimates of the solutions to these problems. 


6.1 A’priory estimates for the weak solutions to the 
dynamic problem B*(r) 

Everywhere below, for simplicity of the presentation omit the index ¢ where 

it does not cause misinterpretation. In addition, the domains 125 (r(@jt)); 

Ne (F(z, t)), Or Get) and the boundary ['(r(a#,t)) we will denote as 

OQ f(t), Q(t), QF(t) and P(t). 

Lemma 6.1. Let r € Wr. 


Then in conditions of Theorem 38 for all weak solutions to the dynamic 
problem B(r) hold true estimates (5.2). 


Proof. We put v in the integral identity (4.2) as a test function . 
One has 


| a, |D(x, v)|?dz = -{ (1 — x)./a, D(a, 0) -F dx < 
Q7() 


Q;(r) 
Pa 4 
Vv a, |D(x, v)|"dx + — 
2; (t) ¥ JOs(t) 


where F= y \/a,, D(z, v) and v is an arbitrary small positive number. 


F \?dx, (6.1) 


Taking v = Pl and considering the uniform boundedness of F with re- 


spect to © (estimate (2.14)) we will get the required estimate (5.2) for the 
summand ,/a, D(x, 0). 

The estimate for the velocity vector follows from the Poincare inequality 
(part 1, § 116, [36]) for the function 0(ke + cy,t) = u*(y,t) = w*(z,t), 
which is equal to zero on the boundary of the domain OF) in the cube eY 
with edge length e: 


if lady <M? i D(y, w*) Pay, 
Y Y 


then _/ ja*|"dy < om | |D(z, w*)|?dy 
eY eY 
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with the constant M independent of «. 
Thus, 


| |/?dy < &? ue f 
ak (r) a (r) 


Dy, ¥)?|dy = 


ve [an (B-B)dy < MAC) 
QE (r 
f 


what guarantees the estimate for the velocity vector. 
To prove the boundedness of the pressure p(a,t) in the space L2(Q7r) we 
represent the identity (4.2) as 


If Vv: (/a, ((F- F) — p° ) -pdxdt = (f,¢), (6.2) 


f=V- (Vou ((F-F)-—p°l), F=x° Ya, D(z, ), F, Fe La(Qr), 
oot zi] 


fen, H=w, (Or), Ho! =L,((0,T);w, ()) 


and 
Ices = sup Lan (ve (P= F) = 91) + Diop) deat] 
ocd D(x, )|lo,a7 S 
IVa ((F-F) — PD llaan < Mp. (6.3) 


The equality (6.2) and the estimate (6.3) mean that a linear functional 


1: Lo(H-!) >R (6.4) 


is bounded (§ 1, Chapter IV, [20]) and and its norm does not exceed the 
value M,, which is equivalent to the estimate (5.2) for the pressure. 


6.2 A’priory estimates for the weak solutions to the 
diffusion — convection problem B*(r). 


Lemma 6.2. The diffusion — convection problem B*(r) has an unique weak 
solution for which the estimates (5.3) — (5.5) hold true. 


Lemma 6.3. The diffusion — convection problem B*(r) has an unique weak 
solution for which the estimates (5.3) — (5.5) hold true. 
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Proof. 'To prove the lemma, we use Theorem 3 (Leray-Schauder fixed point 
Theorem). Let 9% = {u € Lo(Qpr(r)) : |lell20,o0) < M(Mo)}- 

It is easy to see that this set is closed in the metric of space V2(Qy,r(r)). 

Next we define an operator s’ = W(u), s7 is the solution to the modified 
diffusion — convection problem B(r), where the function ~(c) is replaced by 
the function 7 ~(u). 

Let’s call this linear problem as a problem B’(r, u). 

By construction the solution s7 to the problem B’(r, uv) satisfies an inte- 
gral identity 


Feo? 8" (x, to) E(x, ~ s™ (a, 0) €(a, 0)da+ 


i ee ate: (do Vs" —rdY(u))dxdt = 0, (6.5) 


which holds true for all t),0 < t) < TV and for all smooth functions €, 
vanishing at the boundary S'U S? for t > 0. 

The linear problem B’(r, wu) has an unique weak solution s7 from the space 
V2(Qy,r(r))and estimates (5.3) and (5.4) are fair for its solutions. Due to the 
simplicity of this statement, we were unable to find publications with proof of 
the formulated statement, but its solvability is simple enough to prove using 
the Galerkin method and the methods developed in the monograph [25]. 

Estimates (5.3) and (5.4) for s7 and Theorem 2 show, that the operator 
W’7(u) is completely continuous and maps the set Nt into itself. 

For 7 = 0 the problem B°(r, u) has a solution s°. Therefore, the operator 
W7(u) has at least one fixed point s! = c which will be the weak solution to 
the modified diffusion — convection problem B!(r, c). 

Only now that can we prove the estimate (5.5) for the solutions c to the 
modified diffusion — convection problem B(r). 

We put in the identity (4.6) c=u+c* and € = ut = max{u, 0}. 

By definition ut > 0, utu=|ut/??, Vu-Vut =|Vut|?, ut =0 at the 


+ 
boundary St US? and owe = OM ie 


Taking into account the formulas (4.5) we arrive at a chain of the equal- 
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ities: 


| (c* + u(@, to) ut (a, toda = | (cut (a, to) + |u(a, to) |?) dx = 
Qf (to) 


OQ (to) 


to + 
| i (ute) drdt — +h = 
0 JOF(t) t 


ra a) 1 
| | —(cut += |ut)P\dedt + -—bh= 


1 
| (cut (a, to) + = |u*(w,to)|?)de + Ip +h — be. 
(to) 2 


to to 
I, =do i, Vut-Vudzrdt = do fi i |V ut |? dedt, 
0 YO s(t) 0 YOQ;() 
to to 1 
In= i Vul-ow(ul+c*)drdt, Ip = | ip (cut + = |ut|?) DS, sin y dodt. 
o Joye) o Jrit) 2 


Since Di, siny > 0, then we finally get the inequality 


to 
: : |u*(a, to)|\"dax < -{ Vul-ouv(ut c*)dadt, (6.6) 
2 Jos (to) 0 Jose) 
from which it follows that ut (a,t) = 0 almost everywhere in 0,7. 

In fact, otherwise on the set Q, where u > 0 (or where c > c*) the left 
hand — side of the inequality (6.6) strictly positive, whereas the right hand 
—side of this inequality is zero. That isc < c*. The case c > 0 is considered 
similarly. 

Estimates (5.5) the solutions to the modified diffusion — convection prob- 


0 
lem B*(r) show, that in the identity (4.6) w(s) = s+ 7 That is, this prob- 


lem coincides with the diffusion — convection problem B*(r) and for extended 
functions 0° and c the identity 


[ xo, 2)e(.to) x(r(esta), Zar — f x(r(@.0),2) (a) E(x, 0)de-+ 


E 
to _€a Saree 
/ [xeBl- BEL VE. (VE —F (c + =)))dedr =0 (6.7) 


fulfilled for all smooth functions €, vanishing at the boundary S$! U S? of the 
domain Q for t > 0. 


Remark 3. In particulary, the identity (6.7) means that 


ac E Le((0, 7); Wy '(Q)) and 
ae 
1 < c): : 
lle < MOM.) (6.8) 
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7 Proof of Theorem 5. 


7.1 Homogenization: the choice of convergent subse- 
quences. 


Estimates (5.2) — (5.4) allow you to select converging sequences as ¢ > 0 
(for simplicity, we will leave the same indexes) 


vv, #'S V, (1—°)0* converges strongly in L2(Qr) to zero; 


ex’ D(z, 0°)" D(y,V), Be +p, pe’ P; 
Oo € 
Ja, F converges strongly in L2(Q7) to zero; 


O° 


U 


© converges strongly in L2(Q,) to zero; 


c*converges strongly in Le(Qr) to c, y* c* ps MC; 
Ver Veoh Ve, (24) 


where functions V = V(a,t,y), P = P(a,t,y) and C = C(a,t,y) are 1 —- 
periodic in the variable y, 


v, pELa(Qr), ce Wy°(Qr), P, V,C, D(y,V), Vy C € Lal, Qr x Y) 
and 
|v ll2,.07 + |lplloor < Mp, 
|Plloyxar + ||Vloyvxa- + |l D(y, V) Iloyxor < M;,; 


0 < c(a,t) < c*, (a,t) € Or, 
IV clln0r + |V Cllayxar < Me. (7.2) 


7.2 Homogenization of the problem B*(r) 


Lemma 7.1. In conditions of the Theorem 4 the macroscopic continuity 
equation 


V-v=6— (7.3) 


and the first boundary condition in (4.16) for the velocity vector hold true in 
the domain Qr. 
In (7.3) 
4 
mr) = i x(x, t,y)dy =1-—57r°. 
y 3 
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Proof. The limit as « — 0 in the identity (4.3) with smooth test functions 
w, vanishing at the boundary S'US? of the domain 2 for t > 0, gives us the 
integral identity 


= oy J ee 
o= | f (-omE+v-vuydear= f f (6 Ot w+vu-Vw) dxdt. (7.4) 


The identity (7.4) means that differential equation (7.3) and the first bound- 
ary condition in (4.16) are satisfied in a usual sense. 


Lemma 7.2. In conditions of the Theorem 4 the limiting functions v and p 
are bounded in the spaces L2(Qr) and W3°(Qr) correspondingly: 


Ile ]2.0n + IV pllaor < Mp; (7.5) 


and satisfy the Darcy ’s law of filtration 
1 
vas evn?) (7.6) 
1 


in the domain Qr and the second boundary condition (4.16). 
A symmetric strictly positive definite matriz C’(r) is given by the formu- 
las (4.19). 


Proof. First of all we prove that 


P(x,t,y) = p(w, t) x(a, t, y). (7.7) 
x 
To do that we put in the identity (4.2) yp = € yo(a,t) y,(—), where yo(z, t) 
is an arbitrary smooth function in the domain (27, vanishing at the lateral 
boundary of the domain Qr, and y,(y) arbitrary smooth function in Y. 
Passing to the limit as « > 0 we get 


as ea(e.t( | Ple.ty) Vy: Pi(y)dy) dxdt = 0, 


which is equivalent to equality (7.7). 
The continuity equation and the boundary condition for the velocity vec- 
tor V: 
Vy-V=0, ye Ye(r), V =0,7 |yl=r, O<t<T (7.8) 


are derived in [15]. 
To derive continuity equation we take as test functions w in the integral 
x 
identity (4.3) functions of the type w = € Wo(a,t) vil) and pass to the limit 
as € — 0. 
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To derive boundary condition in (7.8) we consider a two — scale limit in the 
equality (1—x*) #° = 0, which leads to the relation (1—-x(a, t, y))V (a, ty) = 
0. 

Next we use the smoothness of the function V: V € W,°(Y), which 
dictates the boundary condition in (7.8). 

The limit as « — 0 in (4.2) gives us the Homogenized dynamics 
equation 


[fe (/ (41 Diy, V) : Diy, ) —p) V- p dy) deat = 


i 
-f [ ve? eat. (7.9) 
0 .4G 
a 


Now we put in (7.9) p = y(a, t) i(— ), where solenoidal functions ¢,(y) 


are equal to zero at the boundary y(r) of the liquid component Y;(r) of the 
domain Y, y, € W3(Y;(r)), supp y; € Y;(r), v:(x,t) = 0 at the boundary 
S° of the domain © for t > 0, supp y; C Yp(r) and (~;)y, = e:, where 
{e1, €2, e3} is orthonormal basis in R?. Such a choice is possible due to the 


Lemma 2.1. 
ik 
I / (Yo a; — p O90) devat = -{ i. V p’ dxdt, (7.10) 
QA On, 0 JO 


One has 
ay = | Diy, V) : Dy, Prs)dy € La(Qr). 
Yf 


Therefore, a; = op and 


Ox; 


3 
IV plo, < MYO | > a D(y, V) : Dy, Pi ,)dylla0r < Mp, 
i=l Yy 


which is proves the estimate (7.5) for the pressure. 

Since in the identity (4.2) functions yoare equal to zero only at the part 
of the boundary S° of the boundary S for t > 0, then p is obliged be zero 
at the part S' US? of the boundary S for t > 0, that is proves the second 
boundary condition in for the pressure in (4.16). 

Since the pressure belongs to the space W,’(Q,r), then equalities (7.8) 
might be rewritten in the form of the Stokes equations 


ti AyV —-Vyl=Vi(p+p), Vy. V=0, ye Yy; 
Viz,t,y) =9, ye 77). 
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The solution of the last system is given by an explicit formula 


V=-=J(Vi@e)- Vip +p, 


PA es 
where functions V‘(y) are defined by formulas (4.19). 
Thus, 
be 1 
v=(V)y, = = @ei))y,- Vptp’) = Aer (p+ p°), 
i=1 


which is complete the prove of the Lemma. 


Consequence 7.1. The limiting functions u(x,t) and p(a,t) are defined 
uniquely. 


The proof of the last statement follows from the estimate (7.2) for the 
difference of two possible solutions to the system of Darcy’s system of differ- 
ential equations (7.3), (7.6), completed with boundary conditions (4.16). 


7.3 Homogenization of the diffusion — convection prob- 
lem B*(r). 


Lemma 7.3. [n conditions of the Theorem 4 the diffusion — convection prob- 
lem H(r) has a unique weak solution c € Vo(Qr) such that 


Iclloa- + ||V cllaa, < Me, O< c(w#,t) < c, (a, t) E Or Grell) 


and the integral identity 
| mc(a, to) €(a, to)dx — | mc?(a) E(x, 0)dx+ 
Q Q 


2 ( dg ' 0 
— mc 4+ VE: (do C(r)Ve— vu (c+ =))) deat = (,. (72) 
o Ja ot 6 
is fulfilled for all smooth functions € vanishing at the boundary S* US? of the 
domain OQfor t > 0. 
Symmetric strictly positive definite matrix C°(r) is defined by formula 
(4.19). 


To prove identity (7.12) it is enough to pass to the limit as « > 0 in the 
identity (4.4). 

Estimates (7.11) follow from estimates (5.3) and (5.5). 

Properties of the matrices C’ and C* have been studied in [15] for the 
case of the fixed structure of the pore space. 
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In our situation these matrices will only be positively definite for all 
r(a,t) > 0. That is, the result we need is a uniform evaluation from below 
of the eigenvalues of the matrices C’ C°, which is valid for all r(a,t) > 0. 
Therefore, we have to prove additionally the desired statement. 


Lemma 7.4. In conditions of the Theorem 4 symmetric strictly positive defi- 
nite matrices C’(r) and C°(r) are infinitely smooth with respect to parameter 
re 


(C°(r)-¢) -¢ + (Cr) -¢)-¢ & Hm >0 (7.13) 
with strictly positive constant v9 for all € € R°, |¢| = 1. 


Proof. First of all we will consider the matrix C’(r), prove the solvability of 
the problem (4.19) for r(a,t) = 0 and study the properties of its solutions 
V'(r;y) for all r(a,t) > 0. 

Recall,that matrix C’(r) is defined by solutions V‘(r; y) to the boundary 
— value problem (4.19) 


A,V' —V,II' = -e’, V,-V' =0, ly| >, 
Viy)=0. Ilan [| Wedy=o (7.14) 
Yy 
as 
3 4 * 
C'(r)=25 (Vy, Be. 
i=1 


It is easy to see that elements of the matrix C’(r) are infinitely smooth with 
respect to parameter r. 

Multiplication of the dynamic Stokes equation in (7.14) by V? and inte- 
gration by parts over domain Y;(r) give us 


(D(y, V’) : Diy, V')y, = (VW? -e’)y,, i, = 1,2,3. (7.15) 
Relations (7.15) show, that the matrix C’(r) is symmetric. 


These equalities for 7 = 7 and the Poincare - Friedrichs inequality (§ 116, 
chapter III, Vol. 4, [86]) imply the estimate 


i D(y, V') Pay = |( 
Yy(r) 


1 1 
V | |VPdy + — Mr & v M3 | ID(y, V") (?dy + me) 
¥j(r) aH ¥4(r) au 


Vidy) -e'| < 
Y¢(r) 


for an arbitrary positive v with some constant M3 in the Poincare - Friedrichs 
inequality independent of r. 
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1 
Taking v M3 = 5 we obtain 


ie (\V'? + |[D(y, V’)|?)dy < M3m(r). (7.16) 


Functions V’ and D(y, Vv") are continuous with respect to the parameter r 
lim V =V5, lim D(y, V’) = D(y, V4), 


where 
AyVo — Vyll, = —e:, Vy: Vo =9, ly] > 0, 
Vi(r:0) =0, | i (r;y)dy = 0 (7.17) 
Y 


Let now consider arbitrary constant vectors ¢ = (C1, G2, ¢3) and 7 = (1, "2, 73) 
from R® with unit norms and put 


=DGV"(y), nly) = 2 mV"(y) 


fm) = (Cr) -n) «n= [ [Dr zalPav 


The continuous function f(7) reaches its minimum f°? = f(7°) on the sphere 
In| = 1. 

Suppose that f° = 0. Then, D(y, 2,0) = 0, which which is possible only 
if Z, is a linear function of the variable y. 

On the other hand — side the boundary — value problem (7.17) implies 


A,z — Vyll =—n°, Vy 20 =9, lyl > 9, 
z,,0(0) = 0, i IL,o(y)dy = 0. 
Y 


Therefore II,o = 7° - y everywhere in the cube Y. Since the function II,0 is 
periodic in y, it is possible only if 7° = 0. 

The obtained contradiction proves our statement. 

Let us now derive formulas (4.19) for the functions C’(y). 

First we consider the limiting identity, resulting from identity (4.4) as 
é — 0 with test functions € = €(a,t): 


0. 6 
o- f i (—meS +a VE- (Vet V, Cdy — v(c+ =) )dxdt = 


Y;(r) 


[Lc med + dy VE (Ver f ” Cdy + f))dxdt = 
Or 


i: A ~meS S 4 dy VE: as. Ve+f))dedt. (7.18) 
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Here 


do F(a, t) = (0 +5) = (fi fa fs) 


Cr) = +f VyCdy =14+ Co(r). 
Yy(r) 
Next as test functions we choose functions € = € (x,t) &(=). 
The limit as ¢ — 0 gives us identity 


[fel Vyei- (Vy C + f)dy) dxdt = 0, 
0 Ja Yy 


where from 


OC 


with n = z 
r 
As usual, we use the representation 


3 


C(r;a,t,y) = >) C'(r:y) fila, t). 


i=l 
One has 
‘i i ac" i 
Vy (VyC’ +e) =0, yey, ae ea, ly|=r. (7.20) 


Multiplying the differential equation for C* in (7.20) by C? and integrating 
by parts over domain Y;(r) we arrive at the equalities 


VyC'- Vy Co dy + VC? -e'dy=0, 7,9,=1,2,38. (7.21) 
Y;(r) Y¢(r) 


In particular, 
i IV cPay+ f Vie re dy =U. 7.= 1,253, 
Yy(r) Yp(r) 


whence follows the a’priori estimate 


3 3 
>| Vcd =|) f V,Ci-e'dy| < 
a CE) j= 2 Ye(r) 


3 
VC’? dy + = m(r); 
fai I ¥4(r) 2 


3 
|, IV, C2. dy < Mum(r), 
j=1 2 Ye (r) 
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with the constant M, independent of the choice of the function r € Nr, and 
unambiguous solvability of the boundary value problem (7.20) in the space 
W3(Y(r)) for all r > 0. 

Besides, from (7.20) follow equalities 


(Co(r) +e) -e7 + Vey, C dH 0, 7 = 182.3, (7.22) 

Yz(r) 
That mean the symmetry of the matrix C§(r) and hence the matrix C‘(r). 
Strict positive definiteness of the matrix C°(r) is proved completely anal- 
ogous to the strict positive definiteness of the matrix C’(r), which completes 
the proof of the Lemma. 


Consequence 7.2. We can always assume that the matrices C’(r) C*(r) 
are diagonal. 


To prove this fact, it is enough to refer to the Theorem 5.28 (§ 6, Chapter 
5, [38])on the reduction of two symmetric matrices, one of which is positive 
definite, to a diagonal form by means of an orthogonal transformation of the 
coordinate system. 

In this case, a positive definite matrix is reduced to a unit matrix. 

Everywhere below we will assume that our Cartesian coordinate system 
{€1, €2, €3} is an orthonormal system of eigenvectors of the matrix C’(r) = 

3 


and C%(r) = S- c'(r) e; @ e;. Moreover, due to the symmetry of the set Y+ 
i=1 
regarding orthogonal transformations 


1 
en=Scn=e 0) = Ss 
do 
Consequence 7.3. The solution {v, p, c} to the problem H(r) satisfies the 


following initial boundary — value problem 


(r) > sq = const > 0. 


1 O 
v=-—V(ptp), V-v=d =, (a,t) € Or, 
1 Ot 
Op _ 0 = 1 2. 
apo (ee) eo, pO. (ete SUS. 17-23) 


—(mc) =V-(s(r)-Ve—v(et+ “)). (x,t) € Or, 


c= 0°, te 5, US, 
Oc 
= 
where 7 is the unit normal vector to the boundary S%. 
The initial boundary — value problem for the limiting concentration c 
is understood in the sense of distributions as the integral identity (4.17), 
completed with the boundary condition (7.24) at the boundary $7. 


0, (w,t) € S%, c(w,0) =c°(x), 2 EO, (7.24) 


AT 


7.4 Differential properties of solutions to the problem 
H(r). 


Lemma 7.5. In conditions of the Theorem 4 the solutions v and Vp to the 
problem (7.23) are from H™? (Qr) and the following a’priori estimate 


Vol? + [vl < M(M,) (7.25) 
hold true 


Proof. The function p(x,t) satisfies in the domain 2 for t > 0 the following 
boundary — value problem 


Om 
Ap= = — p56 — 
P= od = va 
p(a,t)=0, xE€S'US’®, t>0, oP (a@,t) =0, gees, t>0; 
f(.,t)€H = [0,7]. (7.26) 


Let p(a,t) be the extension of p(a,t) across the boundary S = OC into the 
domain Q = {a € R®: |zx,| < 1, i = 1,2,3} in the direction of the normal 
vector to S in an even or odd way. Function p in each unit cube from Q \ Q 
coincides with the shift of the function +p to this cube. 

At the same time, the extension p(x, t) satisfies in the domain Q for t > 0 
the Poisson equation 


Ap=f. (7.27) 
At the boundary $$ = {x2 = +1} U {x3 = +1} of the domain Q 
Op 
a = 2 
ae 0 (7.28) 
and at the boundary Sj = {a, = (+1)*} Q 
p= 0. (7.29) 


Due to the boundary condition (7.29) and the inequality 


max || f(.,t)|laq < M(M,), 


0<t<T 


for the solutions to the boundary — value problem (7.27) — (7.29) holds true 
the inner estimate (Lemma 9.12 §11, Chapter 9, [37]) 


max ||VA(.,t)[[0 < max ||f(.Dllae < M(q,Mp), (7.30) 


0<t<T 0<t<T 


where q > 1 is any. 
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Let y= i Then 


IVAlnn + (VP Ror S max IVa < MO, M,) (7.31) 


[25] (Theorem 2.1, §2, Chapter II), which proves the boundedness of the 
pressure gradient and its Helder norm with respect to the spatial variables. 

To estimate the Helder norm of the pressure gradient with respect to time 
time, we consider the differences 


D,B(a, t) = ~ (D(a, t + h) — B(x, t)), 


Deir ~ (Fla,t +h) — Flw,t)). 


These finite differences satisfy the Poisson equation (7.27)) in the domain Q 
and boundary conditions (7.28)) and (7.29) at the boundary Q. 
Therefore, for D,p the estimates (7.30) and (7.31) are fair: 


IDiVB(@,t)| < |DnVPlo, < M(a, M,)- 
The limit as h + 0 in the last inequality proves the estimate 
(3) 
(Vp)io7 < M(q, Mp), 


which is together with the estimate (7.31) and Darcy’s law (7.23) complete 
the proof of the Lemma. 


Lemma 7.6. In conditions of Theorem 4 the limiting function c(a, t) belongs 
we 2 = 
to the space H2+Y =" (Gp) and the estimates 


et? < M(M.) (7.32) 
and 
0<c(a,t) < c*, (x,t) EQ (7.33) 


are fair for her. 


The estimate (7.33) follows from the estimate (5.5) after the limit as 
ee al 

The estimate (7.32) follows from the estimate (7.33) and the local esti- 
mates for classical solutions to the equation (4.18) [25] (Theorem 10.1, §10, 
Chapter IV). 

The existence of the classical solutions to the diffusion — convection prob- 
lem H(r) follows from the estimate (7.32) and Theorem 5.2 5.3 [25] (85, 
Chapter IV). 
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7.5 Properties of operators F’(r), F?(r) and F“(r). 


By construction, operators F°’(r), F?(r) and F° n ) act from the set SW into 
the spaces H”?(Q7), H”?(Q7) and H?+” aa(9) QQ) accordingly. 
Takes place the following 


Lemma 7.7. In conditions of the Theorem 4 the operators F’(r), F?(r) and 
F°(r) are Lipschitz continuous 


IF’ (r1) — F°(r) | + [FP (r1) — FP(r2) |G + 
IF°(r1) — F°(r2)|G? << M(Mo)|ni — ral SF”. (7.34) 


The proof of this statement is standard and expresses the well-known 
fact of continuous dependence of solutions to linear elliptic and parabolic 
equations on coefficients. 


8 Proof of 5. 


It is easy to see, that the operator F(r), defined by formula (4.27) satisfies 
the Lipschitz condition. Moreover, for some small time interval (0,7) it is 
compressive and displays the set tr into itself. 

t 


In fact, let R;(a,t) = F(r;) = | c(x,T)dr, 1 = 1,2, where c; = F°(rj). 
0 
One has: 


F(r)|G? < T,M(M), 
Fn) — Er Hla » < T, M(Mo) |r1 — rele”. 


0 < F(r)\(z,t) < Tc’, 


That is, on the interval (0,7;), where 


M(Mo) 
2 


1 
(ie mint } 
the operator F(r) is compressive and displays the set Mtv, into itself. 

Banach ’s Theorem (Theorem 1, §4,Chapter I, [20]) guarantees us the 
existence of the unique fixed point r*(a,t) from the set My, and, that is, the 
validity of the Theorem 5 on the time interval (0, 7;). 

Next we put r;(x) = r*(a,7,) and consider the problem B*(r) on the 
interval (7,7), where instead of function ro(a) we will take the function 
r,(a), and instead of function 7(a#,t) we will take the function 7)(a,t) = 
max{0, ri(a) — r(ax,t)}. 
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Let also Qn, 7) = Q x (L1,T) and 


coe 2y 1 
Merry = {re H+" (Quer, 1), Te, = 0.0 ea 5° 
Or, 2 
—@ < ap (et) < 0, rile < Mo}. 


After that, as before, we will define solutions {v1, p1, ci} and {rj, v1, pi, a} 
of the problems H(r) and H accordingly , at intervals (71,7) (71,7). 
Repeating the process we will find the time intervals (7;,7k41), k = 
1,2,3,... and a function r*, is equal to the function rf on the intervals 
(Th, Te11), — 1, 2: baad 
Obviously, these functions solve the problem H on the intervals (0, 7;). 
If im Tr = co, then the theorem is proved. 


If jim T, = T* < «© and r*(x, T*) is nonzero on some open set of , then 


by virtue, of the obtained estimates of the solutions of the problem H we can 
calculate the limits of solutions as t + 7* and then will find the solution of 
the problem H on the interval (7*,7* + 6) with some small 6 > 0, which 
contradicts our assumption. 

Thus, the process will be terminated only if r*(a,7*) = 0 on the set Q. 
For t > 7™ the liquid will completely fill the domain 2 and its motion will be 
described by the Stokes equations, which completes the proof of the theorem. 
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